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Abstract. Jet quenching measurements using leading particles and their cor- 
^ ' relations suffer from known biases, which can be removed via direct recon- 

1} . struction of jets in central heavy ion collisions. In this talk, we discuss several 

modern jet reconstruction algorithms and background subtraction techniques 
that are appropriate to heavy ion collisions. 
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1. Introduction 



In high energy collisions, during the hadronization process, colored partons from the 
hard scatter via soft quark and gluon radiation form a "spray" of colorless hadrons 
' also known as "jets" . Jets are expected to reflect kinematics and topology of partons 

and they are the experimental signatures of quarks and gluons [HI [21 [31 H] ■ The final 
states of hadrons can be grouped into jets via various algorithms. All the algorithms 
that are appropriate for hadronic collisions can be used to crosscheck and estimate 
the systematics of the measurement of jets and the fragmentation functions. At 
Tevatron, the inclusive jet cross section is measured for 20 orders of magnitude via 
sequential and cone algorithms and these measurements are consistent with the NLO 
pQCD prediction from CTEQ 6.1 [[51 [5]. The first direct measurement of inclusive 
mid-rapidity jets at RHIC are performed in polarized p + p collisions at y/s = 200 
GeV with STAR experiment [[7|. For this measurement, the neutral particles are 
measured with the barrel calorimeter of STAR with full azimuthal (0 < <j) < ^ 7r ) an d 
partial rapidity (0.2 < 77 < 0.8) coverage. The neutral particles together with the 
charged particles that are identified with the time projection chamber (0 < <f> < 2n, 
\rj\ < 1.3) are combined into jets with a mid-point jet cone algorithm where split- 
merge step is required. This inclusive jet cross section measurement from RHIC at 
\/s = 200 GeV polarized p + p collisions also agrees well with the NLO pQCD. 
Jets can be used as direct probes of partonic phase in the heavy ion collisions at 
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the RHIC and soon at the LHC [ 8 , 9J . Modified leading logarithmic approximation 
(MLLA), which provides a good description of vacuum fragmentation, introduces 
medium effects at parton splitting [ [TO] . In this approximation, fragmentation of 
jets is strongly modified for hadrons with transverse momentum of px = 1 — 5 GeV. 
The study of fragmentation of jets to probe the high density medium created at 
heavy ion collisions is similar to the study of deep inelastic scattering experiments 
as a probe of spatial structure and the forces introduced. During the last 8 years 
of RHIC operations, instead of direct reconstruction of jets due to the challenges 
in understanding the underlying high multiplicity events, measurements of high 
Pt hadron suppression are studied via di-hadron correlations and nuclear modifi- 
cation factors (Raa) [[HI IS]- While these partonic energy loss measurements are 
consistent with pQCD based energy loss through medium induced gluon radiation 
calculations, surface emission bias limits sensitivity to essential q measurement and 
they can only provide a lower bound to the initial color charge density [ [TTJ . Full 
unbiased jet reconstruction (i.e., reconstructing partonic kinematics independent of 
quenched or unquenched fragmentation details) is required in heavy ion collisions. 
Jet shapes, fragmentation functions, and energy flow will increase our sensitivity 
to quenching. In the following sections, we discuss several modern jet reconstruc- 
tion algorithms that are appropriate for heavy ion collisions and their background 
subtraction techniques. 

2. Algorithms Used 

See [ I12| I13[ I14j and references therein for an exclusive overview of many tra- 
ditional jet algorithms used in high energy collisions. Here we discuss only two 
algorithms; seeded cone and sequential recombination, and their background sub- 
tractions. As in cleaner e + + e~ or p+p collisions, in heavy ion collisions the chosen 
jet reconstruction method should be theoretically and experimentally consistent. 
For example from the theory side, the algorithm should have low insensitivity to 
hadronization, radiation, splitting and be equally well defined at hadron and parton 
level. The larger LHC luminosities (20 to 200 collisions in a detector) require that 
the traditional jet algorithms have to be improved to resolve pile up of events. The 
underlying pile up event can be subtracted with the clustering algorithm and the 
consideration of the jet area [[HI [16] . These subtraction techniques can also be used 
in a high multiplicity heavy ion environment where the background subtraction is 
required. As Weinberg says "Quark and gluon jets (identified to partons) can be 
compared to detector jets, if jet algorithms respect collinear and infrared safety." [ 
[4] . However, the cone algorithms with seeds can be used to estimate the systematics 
of the background subtraction in heavy ion collisions. The algorithm should also 
be detector independent and should allow the combination of particles detected in 
various detectors such as barrel calorimeter and time projection chamber for neutral 
and charged particles. As the physical energy measurement degrades the resolution 
somewhat with the detectors, the choice of algorithm should be aimed to minimize 
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the additional resolution effects. To handle a large number of events with high 
multiplicities, the computational efficiency is also required. However, with the im- 
proved grid system and continuously increasing technology of computational speed, 
it is possible to work with slower algorithms [ [17] . 

2.1. Cone 

The main principle of the cone algorithm is to combine particles in rj — <j> space 
with their neighbors within a cone of radius R (R — y A(f) 2 + A77 2 ). Variation of 
the algorithms such as seeds (requirement of minimum energy), splitting, merging, 
iteration and so on exists to optimize the search and effectiveness of jet finding. 
Depending on these specific requirements of the algorithm, infrared and collincar 
safety is not a guarantee. Selection of the cone size and seed can also cause trigger 
biases. However, cone algorithms are very intuitive and are based on the picture 
that a jet consists of a large amount of hadronic energy in a small angular region. 
That is why they have been used as a primary tool to identify jets at hadron colliders 
starting since the early 1980s. We propose to use a simple cone algorithm with a 
seed cut and without the split-merge step to study and estimate the background 
subtraction of the complex heavy ion collisions. 

Figure [l]-a shows a comparison of the total energy within a cone size of 0.7 for 
the two pt cuts for the most central and peripheral Au + Au collisions of simulated 
Hijing events [HI]- For larger px cut (pr > 1 GeV) the total energy in a jet 
cone is reduced significantly (~ 50%) with respect to no pr cut, lowering the total 
background energy collected in a jet cone. Figure [T]-b shows a similar comparison 
for the total energy in a jet cone within two different cone sizes when a fixed px 
cut is applied. Reducing the cone size also reduces the net energy collection in the 
jet cone as can be seen as a shift towards the left in the energy distributions. By 
using a smaller cone size and a larger pr cut, it is possible to suppress the heavy ion 
background. The smaller cone size also implies a reduction of the total measured 
jet energy. However a small cone size of R = 0.4 can be appropriate, as ~ 80% of 
the jet energy is observed to be within R~ 0.3 for 50 GeV jets in the Tevatron [ 
6 . This assumption is true only if the broadening of the jet fragmentation due to 
quenching in the medium is small. Even with the larger px threshold, additional 
background subtraction is needed. One possible way is to subtract the rest of 
the background from the averaged out-of-jet-cone event-by-event for the selected 
cone size. Quantitative studies are needed when selecting regions of subtraction 
due to r\ and <f> acceptance, R and px cuts, fluctuations and elliptic flow observed 
in real heavy ion collisions. The overall energy resolution, detector efficiency and 
acceptance can be estimated and corrected by simulated quenched or unquenched 
fragmentation of jets that are embedded in real heavy ion events. 

To study the effect of the background in jet reconstruction for the LHC colli- 
sions, quenched jets are generated and reconstructed with a simple cone algorithm 
without any iterations or split-merge steps in p + p and Pb + Pb collisions [ [TH1 HO] . 
In these studies, the Pb + Pb background is observed to have a minor effect on 
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Fig. 1. Comparison of the total energy (a) for the two pr cuts (0.1 and 1 GeV) 
for a fixed R = 0.7 and (b) with two cone sizes (R = 0.4 and R = 0.7) for a fixed 
Pt cut. For these plots Hijing events are generated and separated into most central 
and peripheral Au+Au collisions at ^snn = 200 GeV. 



the reconstructed jet-energy, ft is also observed that the heavy-ion background is 
manageable as the differential cross section of reconstructed jet energy with and 
without Pb + Pb background seems to agree well for jets with Et > 50 GeV. 

2.2. Sequential Recombination 

The sequential recombination algorithms, also known as the kx algorithms, have 
been used extensively in the Tevatron as they are collinear and infrared safe [ 
[21] . In these type of clustering algorithms, arbitrary shaped jets are allowed to 
follow the energy flow resulting in less bias on the reconstructed jet shape than 
with the seed base cone algorithms [ [TB] . The Ut algorithms include various steps 
such as clustering objects close in relative pr, merging these particles into a new 
cluster when the momentum satisfies a predefined minimum, and repeating steps 
until all clusters become jets. Details of this type of algorithm can be found in 
[ [6] [21] and references therein. While the steps in combination of particles in 
sequential recombination algorithms might not be as intuitive as the ones from the 
cone algorithms, the framework of clustering algorithms such as fastjet are available 
and easy to implement on experimental data [ [23]. With the assumption that 
underlying event and pile-up are distributed uniformly in 77 and </>, an active area 
(Aj) of each jet in the fastjet is estimated by filling an event with many very soft 
particles and then counting how many are clustered into a given jet. Study of Pt /Aj 
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of jets determines the noise density and can be subtracted from the measured jet 
energy on an event-by-event basis to correct for the real jet energy. This correction 
is observed to recover most of the simulated jet's energy when they are reconstructed 
in pile up and heavy ion backgrounds [ [16] . A new reverse clustering type called 
anti-fc-r is also available within the fastjet algorithm framework [124). The anti-fey 
algorithm provides a fast infrared and collinear safe replacement of the traditional 
cone algorithm and fixes the irregularities in the boundaries of the final jets due 
to the soft radiation observed in regular kx type algorithms. Since jet areas can 
be subtracted as in the kx case, this reverse sequential clustering algorithm is also 
appropriate for the heavy ion environment and is essential to implement on heavy 
ion events as a tool to study the systematics of jet reconstruction. 

3. Conclusions 

The LHC is expected to start at the end of the summer 2008 and will be a jet factory. 
At these high energies, the heavy ion background will be larger but also there will 
be copious production of jets that will be sufficient enough to get above the heavy 
ion background [|25]. High cross-section rates of jet production provide sufficient 
energy lever-arm to map out the energy dependence of jet quenching. At RHIC the 
backgrounds are smaller but also cross-sections of high momentum jets are smaller. 
Jet structure changes due to energy loss and the additional radiated gluons are 
expected to be observed at LHC with respect to RHIC. It is essential to reconstruct 
jets at LHC and RHIC to measure the energy dependence of the medium induced 
energy loss to map out all energy levels. The detector set-ups like calorimeters 
and time projection chambers are available from both the STAR experiment at the 
RHIC and the ALICE at the LHC. The clustering and cone algorithms provide 
the necessary systematic study of infrared and collinear safe jet reconstruction and 
of the variation of techniques applied in subtraction of large heavy ion underlying 
events and fluctuations. While the energy resolution depends on fragmentation 
functions of hard scattering cross sections, the question of quenching can only be 
addressed by theoretical input. 

Appendix 

First jet reconstruction measurement in most central Au + Au collisions at y / sjv r /v = 
200 GeV from Year 2007 run at the RHIC and the systematic study of utilizing both 
cone and cluster algorithms are presented during the Hard Probes meeting in June 
2008 at Ilia da Toxa Galicia in Spain [I2B1 [27]. 
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